The Near— Threshold Production of </> Mesons in pp Collisions 
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The pp~*pp(f> reaction has been studied at the Cooler Synchrotron COSY- Julich, using the in- 
ternal beam and ANKE facility. Total cross sections have been determined at three excess energies 
e near the production threshold. The differential cross section closest to threshold at e = 18.5 MeV 
exhibits a clear 5*-wave dominance as well as a noticeable effect due to the proton-proton final state 
interaction. Taken together with data for ppw-production, a significant enhancement of the <j>/uj 
ratio of a factor 8 is found compared to predictions based on the Okubo-Zweig-Iizuka rule. 



PACS numbers: 25.10.+S, 13.75.-n 

Meson production near threshold has the potential to 
clarify important questions of hadron physics in the non- 
perturbative regime of quantum chromodynamics due to 
its comparatively simple scheme of interpretation. The 
production of light vector mesons, p(770), w(782) and 
0(1020), quark anti-quark states with their spins aligned 
(J p =l~) and without open strangeness, has been inves- 
tigated with both hadronic and electromagnetic probes 
in order to study production mechanisms |l|, coupling 
constants 1], modifications in nuclear medium |2| and 
in particular the so called Okubo Zweig Iizuka (OZI) 
rule |3j . This rule states that processes with disconnected 
quark lines between initial and final states are suppressed 
compared to those where the incident quarks continue 
through to the exit channel. As a result, the production 
of ideally mixed </>-mesons (quark content ss) in a reac- 
tion A + B — ► 4>X is reduced compared to A + B — > ujX 
[ui is a linear combination of uu + dd) under similar kine- 
matical conditions. Taking into account deviations from 
ideal mixing between singlet and octet vector mesons, 
Lipkin predicted a ratio of single <j> to lo production of 
R(j}/ui = 4.2 x 10~ 3 = Rozi llll]- However, strong en- 
hancements of the experimental R$i w compared to i?ozi 
have been observed (an overview is given in Ref. |||), in 
particular in pp annihilations, where R^i^ can be as large 
as ~ 100xi?ozi 7] ■ Here a strong correlation of the <j> me- 
son yield with the spin-triplet fraction of the initial state 
was found |g], and this in part motivated the suggestion 
of a polarized internal strangeness component in a po- 
larized nucleon |9(. However, alternative exp lanations, 
such as two-step kaon-exchange models [iJJ, UJj, have 



also been advanced. Since vector-meson production in 
close-to-threshold pp — > ppV reactions must proceed via 
the spin-triplet entrance channel, the investigation of the 
cross section ratio a(pp^pp(f))/a(pp^ppuj) at small ex- 
cess energies e should provide a clean way of investigating 
possible violations of the OZI rule. 

Total cross sections for w-production in proton-proton 
collisions have been measured in a range of excess-energy 
e from a few MeV up to several GeV [lJjll^lllJI, whereas 
data for pp<j> are very scarce. Two total cross sections of 
cj) production have been obtained for e ~ (2 — 4) GeV, 
but with rather limited accuracy |l,4 |lfj. At low ex- 
cess energy, a single measurement of total and differen- 
tial cross sections has been made by the DISTO collab- 
oration at e = 83 MeV [17j. In combination with the 
u cross section of COSY-TOF at e = 92 MeV 0, this 
yields R^/^ ~ 7xi?ozi- The differential distributions 
from DISTO indicate that cj> production at that energy 
proceeds dominantly via the 3 Pi {pp) entrance channel, 
though other partial waves do contribute significantly. 
To clarify this, it is crucial to extend the measurements 
to such small excess energies that only the lowest partial 
waves can contribute. Such measurements have become 
feasible at the internal proton beam of the Cooler Syn- 
chrotron COSY at the Research Center Julich, using the 
ANKE target and detector facility. Here we report on 
the results for <f> production in proton-proton collisions 
at three beam momenta, corresponding to excess ener- 
gies of e = 18.5, 34.5 and 75.9 MeV. 

ANKE is a magnetic spectrometer [l^J situated at the 
internal beam of COSY. It comprises three dipole mag- 



nets Dl — D3, which guide the circulating beam through 
a variable chicane. The central C-shaped spectrometer 
dipole D2, with a maximum field strength 1.6 T, is placed 
downstream of the target position. D2 is used to sepa- 
rate the reaction products from the circulating beam, de- 
flecting them towards charged-particle detectors on the 
left/right side of the beam for negative/positive charges. 
The hydrogen cluster-jet targ et used provided areal den- 
sities of ~5xl0 14 cm~ 2 |ljj • The average luminosity 
during the experiment was determined through the si- 
multaneous measurement of pp elastic scattering. By de- 
tecting one fast forward-going proton (i? = 5.0° — 8.5°) 
in appropriate detectors, elastic events were easily sepa- 
rated from background. Taking the corresponding cross 
sections from the SAID database |2JJ], luminosities be- 
tween (1.5 — 3.2) x 10 31 cm _2 s~ 1 were determined. The 
uncertainties in these contribute 4% (at e— 18.5 MeV), 
6% (34.5 MeV), and 9% (75.9 MeV), respectively, to the 
final systematic error in the total cross sections. 

The pp — > ppcj) reaction has been studied by detecting 
the K + K~ decay of </>-mesons in coincidence with one 
of the forward-going protons, requiring that the missing 
mass be consistent with that of the non-observed second 
proton. Particle identification relies on time-of-flight 
(TOF) measurements and the determination of particle 
momenta. In the initial step, positive kaons are selected 
by a procedure described in detail in Ref. J21J , using TOF 
between START and STOP scintillation counters of a 
dedicated K + detection system. Secondly, both the co- 
incident K~ and forward-going proton are selected from 
the time-of-flight differences between the STOP coun- 
ters — in the negative as well as in the forward detector 
system — with respect to the positive STOP counter 
that was hit by the K + . These two TOF selections, as 
well as the selection for the K + , were done inside ±3<r. 
The absolute time calibration of all negative and forward 
STOP counters in conjunction with all of the positive 
STOP counters was performed using the abundant 7r + 7r~ 
and w + p pairs. The final selection of the pp — > ppK + K~ 
reaction was made by a ± 3 a missing-mass cut on the 
non-detected proton. This leads to 400 — 1800 identi- 
fied ppK + K~ events depending on the energy. The es- 
timated background inside the proton cut window is 5 % 
(at e=18.5MeV), 12% (34.5 MeV) and 18% (75.9MeV). 

The left panel of Fig. ^ shows the K + K~ invariant- 
mass distributions in the region around lGeV/c 2 . For 
all three beam energies, a clean 4> peak is observed 
at 1.02GeV/c 2 on top of a smooth background of 
non-resonant kaon-pair production. The right panel 
shows the corresponding pp-^pp<f>-^ppK + K~ differen- 
tial cross sections, i.e. distributions corrected for detector 
acceptances (see below). Contributions from misiden- 
tified particles have been subtracted using data from 
outside the proton peak in the missing-mass distribu- 
tions, adding 3% (18.5 MeV), 7% (34.5 MeV) and 10% 
(75.9 MeV) to the final systematic error. Each spectrum 
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FIG. 1: K + K invariant mass distributions at e = 18.5 MeV 
(top), 34.5 MeV (middle), and 75.9 MeV (bottom). Measured 
raw distributions are shown on the l.h.s. while the corre- 
sponding cross sections der/dM are plotted on the right. Only 
statistical uncertainties are shown. The cross sections contain 
a non-resonant part (dashed line, based on four-body phase 
space) and a 0-meson contribution. The solid line is the sum 
of both and includes effects of the detector resolution. 

has been fit with the sum of two contributions. A uni- 
form distribution, based on four-body (ppK + K~) phase 
space, was used for the background, whereas the </> was 
modelled by its natural line shape, folded with a Gaus- 
sian function (a = 1 MeV/c 2 ) to take into account the 
momentum resolutions of the detectors. 

Estimates of the differential acceptance of ANKE have 
been obtained by means of a multidimensional matrix 
Monte-Carlo method, which allows one to determine the 
acceptance independent of the ejectile distributions at 
the production vertex (see Ref. 17]). In general, seven 
degrees of freedom (dof ) are needed to characterize com- 
pletely the ppK + K~ final state, but in our case, a 3- 
dimensional matrix has been used with the following dof: 
(i) the relative momentum of the two final-state protons 
in the (pp) reference frame, (ii) the polar angle of the 
K + -meson in the rest frame of the K + K~ system, and 
(iii) the K + K~ invariant mass. Significant deviations 
from pure phase space can be expected close to threshold 
for both (i) and (ii), due to the final state interaction be- 
tween the two protons, and also the angular distribution 
of the decay of the <j> mesons (see Fig.|3}. The remaining 
dof are contained in the implicit assumption of isotropic 
angular distributions. These assumptions seem in ret- 
rospect to be justified since the resulting Monte Carlo 
simulations reproduce the measured distributions within 
their statistical uncertainties. Each of the three variables 
(i) — (iii) are subdivided into 10 to 30 bins, producing 
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FIG. 2: Total cross section for <f> production in pp collisions 
as a function of excess energy e from this work (circles) and 
DISTO [17f (square). The error bars include both statistical 
and systematic uncertainties. The dashed line shows a phase 
space calculation normalized to pass through the highest en- 
ergy ANKE point, while the solid line, which includes pp final 
state interaction effects, is a fit to all the ANKE data. 



in total several thousand elements, but for all of them 
the acceptance is non-zero. The acceptance corrections 
contribute 10% (18.5MeV), 14% (34.5 MeV) and 19% 
(75.9 MeV) to the final systematic error. 

Using the number of 0-mesons from the fit, the inte- 
gral luminosity for the measurements, and the efficien- 
cies and acceptances of the ANKE detectors, the total 
(/>-meson production cross section has been deduced for 
the three energies, taking into account the branching ra- 
tio in 4> decay of T K+K ~ /T tot = 0.491 0- The results 
are given in Table [I] and plotted as a function of excess 
energy in Fig. [3 Very good agreement is found with the 
DISTO point at e = 83 MeV. The dashed line in the fig- 
ure displays the energy dependence of phase space. When 
this is normalized to the highest energy ANKE point, it 
misses the two lower points by large factors. The solid 
line includes the effect of the final-state interaction (FSI) 
between the two protons in the 1 Sq state using the Jost- 
function method (see Ref. |22j) and scaled such that it 
fits best all three ANKE cross sections. The much im- 
proved agreement here means that it is crucial to include 
the FSI in any description of the data. 

Before discussing the differential cross sections which 
were measured at the lowest excess energy, it is useful 
to note the following constraints. Close to threshold, the 
two final-state protons must be in the x Sq wave, and the 
<j) in a relative S wave with respect to this pair, so that the 
initial two-proton state is the 3 P\_. This in turn requires 
the alignment of the incident (pp)-spin as well as of the 
final 0-meson spin direction to lie along the beam axis 
(see Ref. [lj| for a more detailed discussion). The polar 
angular distribution of the decay kaons in the 0-meson 
rest frame must then display a sin Q% shape relative 
to the beam direction, as is observed for our e=18.5MeV 
data in Fig. Any additional cos 2 0^" contribution, 
induced by higher partial waves, is not visible. In the 
lower part of Fig. |3 we show from left to right the dis- 
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FIG. 3: Differential distributions for e=18.5 MeV. Vertical er- 
ror bars indicate statistical uncertainties and horizontal ones 
bin widths. Upper left panel: da/dfi as a function of the 
cosine of the polar K + in the reference frame of the 0-meson 
relative to the beam direction. The full line is the expected 
sin 2 © shape. Lower panels: da/dQ vs. cosine of the polar 
angle of the 4> meson in the overall cm. system (left), polar 
angle of the emitted proton (middle), and proton polar angle 
relative to the <fi direction (right), the two latter being in the 
(pp) reference frame. Upper right: Dependence of the cross 
section on the pp relative momentum. The dotted line reflects 
pure phase space whereas the solid includes also the pp FSI. 

tributions in (i) the polar angle of the </>-meson in the 
overall cm. system, (ii) the polar angle of the emitted 
protons relative to the beam, and (in) the proton polar 
angle relative to the </> direction. Both proton angles are 
measured in the (pp) reference frame. All three distri- 
butions are consistent with isotropy, as expected for a 
3 Pi-^ 1 5o transition. Finally, in Fig. [3] the differential 
cross section is plotted as a function of the proton mo- 
mentum in the (pp) rest frame. While the phase-space 
calculation (dashed line) misses the data, inclusion of FSI 
for the two protons in the 1 So-state reproduces the ex- 
perimental results. Thus, a clear and significant pp final 
state interaction is observed at e = 18.5 MeV. 

While at e = 83 MeV DISTO also observed the dom- 
inance of the 3 Pi^ 1 Sq transition 17], they did not 
see any indication of pp FSI in their proton momen- 
tum spectrum, and this is consistent with our findings at 
e = 75.9 MeV. Taking both results together, it is tempting 
to ask for the mechanism that suppresses final state in- 
teractions at moderate excess energies: contributions of 
higher partial waves are an obvious conventional cause, 
but more exotic explanations, like a </>A?-resonance (see 
Ref. 22] ), have also been advanced. 

Turning now to the <p/uj ratio, we also present in Ta- 
blc^data on pp — » ppio total cross sections obtained in the 
e range of our measurement |12lll3l ]. The last column lists 



TABLE I: Total production cross section for pp^>-pp4> at our three excess energies (col. 1) compared to pp—yppuj data (12 . Il3f 
(col. 2) at similar excess energies. The last column contains the ratio of cf> to ui cross sections for each line. In all cases the first 
error is statistical and the second systematic. 



4> production (ANKE) 
e^ [MeV] a (tot) [nb] 


u production 
eu, [MeV] a u (tot) [jub] 


<f>/u production ratio 
R Hu x 10" 2 


18.5 43.2 ±2.2 ±5.1 
34.5 84.9 ±6.9 ±13.6 
75.9 188.0 ±19.1 ±41.4 


19.6 ±0.9 1.51 ±0.23 ±0.18 
30.0 ±0.9 1.77 ±0.48 ±0.23 
92 7.5 ±1.9 ±1.5 


2.9 ±0.5 ±0.5 
4.8 ±1.4 ±0.9 
2.5 ±0.7 ±0.7 



the ratios as obtained line-by-line, i.e. at slightly differ- 
ent values of e. Within the stated uncertainties the ratios 
are equal, and we have therefore calculated a weighted 
mean by first fitting and interpolating the u> results to 
our excess energies. This gives 
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R Hu = (3.3 ±0.6) x 10" 
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as compared to an uncorrected weighted mean of the last 
column of Table I, which is about 10 % smaller. Taking 
into account the effects of the finite meson widths on the 
phase space [l2j, |22J changes Rozi by at most 5% at the 
lowest excess energies. 

The production ratio obtained from high energy ppV 
data is ~ (1 - 2.4) x Rozi El El El- Together with 
our findings, this means that there must be a significant 
energy dependence of the OZI enhancement factor [22J . 
which requires more theoretical work to understand its 
origin. In this context let us mention that the exper- 
imental ratio Rh,/ ,,, d educed from itN interaction gives 
(3.2 ± 0.8)xi?ozi |23J, which can be explained in terms 
of the established OZI violation in the <fipir and wpir cou- 
pling [Sj, |24j . The present ratio from near-threshold <fi 
and lu production in pp collisions exceeds this value by 
more than a factor two. It may be a signal for additional, 
and as yet non-understood, dynamical effects related to 
the role of strangeness in few-nucleon systems. 

In summary, we have measured cross sections for (j) 
production in pp interactions at three excess energies, all 
of which are much closer to threshold than previous data. 
The lowest energy result demonstrates the dominance of 
the transition from the 3 Pi (pp)-entrance channel to the 
^So (pp) final-state. Both the total cross section and 
the proton-momentum spectrum indicate a significant 
pp final state interaction. Using data for u-production 
from literature, it is found that R^i^ is about 8x_Rqzi- 
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